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ABSTRACT 

A bright flare from a galactic nucleus followed at late times by a decay in luminosity is often 

considered the signature of the complete tidal disruption of a star by a massive black hole. The flare 
and power-law decay are produced when the stream of bound debris returns to the black hole, self¬ 
intersects, and eventually forms an accretion disk or torus. In the canonical scenario of a solar-type 
star disrupted by a 10® Mq black hole, the time between the disruption of the star and the formation 
of the accretion torus could be years. We present fully general relativistic simulations of a new class 
of tidal disruption events involving ultra-close encounters of solar-type stars with intermediate mass 
black holes. In these encounters, a thick disk forms promptly after disruption, on timescales of hours. 
After a brief initial flare, the accretion rate remains steady and highly super-Eddington for a few days 
at ~ lO^Moyr-h 

Subject headings: black holes, accretion 


1. INTRODUCTION 


Observations of tidal disruption events (TDEs) have 
the potential to unveil super-massive black holes (BHs) 
at the centers of galaxies. In galaxies with quiescent BHs, 
accretion-powered nuclear activity is absent, and TDE 
signatures are, in principle, readily identifiable. The 
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et al.|2010 ,|2013 ). A variety of theoretica 

scenarios have 


from the traditional case of a disrupted main-sequence 
star, to more exotic events involving the total or partial 
disruption of evolved stars such as a white-dwarfs (e.g.. 


Clausen & Eracleous 2011 Krolik fc Piran 2011), red- 


horizontal branch stars ( e.g.,[U 


giants (e.g., [Davies &: King|26 05||Bogdanovic et al. 
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2014), 
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even super-Jupiters (e.g., INikolajuk & Walter||201d ). Of 


particular int erest are the theoret ical models predicting 
ignition (e.g., Rosswog et al.|2009 ) or the ampli fication of 
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number of putative TDE observations grows, more simu¬ 
lations are needed to interpret the wealth of observational 
data. 


T DE mo delin g was pioneered by Rees (1988), Phin- 


ney (1989), and Evans & Kochanek 
les were first in pointing out that. 


These stud- 
:or tfie most likely 


scenario of the disruption of a main-sequence star by a 
10®“^ Mq BH, an UV/X-ray flare followed by a 
decay in luminosity should be expected. The flare and 
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decay are produced by the accretion of bound stellar de¬ 
bris returning to the BH. In particular, the decay 

was singled-out as a ubiquitous signature for the presence 
of massive BHs. Subsequent studies have incorporated 


detailed micro-phvsics (iLodato et al. 12009 

Guillochon & 

Ramirez-Ruiz 2013| 

Rosswog et al. 

2009), considered a 


red-giants dShcherbakov et al. 2013 |Haas et al. 2012 

Kobavashi et al. 1120041 IRosswog et al.i2U09l IBogdanovic 

et al.| 

20141), covered longer dynamical times (Guillochon 

et al. 

2009), and included better descriptions ot gravitv 

tor ul 

tra-close encounters (Laguna et al. 1993 Rosswog 

2010 

Gheng & Bogdanovic||2oi4). 


Coughlin fc Begelman||2014[ [Cannizzo 

1.1 


et al 


lU 
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1990). The complication arises because, in canon- 


ical TDEs, the time between the disruption of the star 
and the formation of the accretion disk amounts to sev¬ 
eral orbital periods of the stellar debris in highly eccen¬ 
tric orbits. This translates into years for a solar-type 
star disrupted by a 10 ®Mq BH. For reference, numerical 
simulations in these cases cover at most a few tens of 
hours. 

Recent papers have addressed the circularizatio n of 
the returning debri s with a variety of methods. 


Sh- 


iokawa et al. (2015) found that the debris circularizes 
at a larger radius than previously thought, and that the 
accumulation of mass in the ensuin g ring is fairly slow. 
Guillochon & Ramirez-Ruiz (2015) considered the self- 
intersection of thin post-disruption streams for an en¬ 
semble of events and concluded that streams typically 
self-intersect at large distances from the BH, leading to 
a long viscous time , hence a long delay b efore the onset 
of rap i d accr etion. Bonnerot et al. (2015) and Hayasaki 
et al. (2015) pointed out the importance of cooling on 
the rate of circularization of the debris, concluding that 
efficient cooling leads to very long circularization time 
scales, hence a delay in the onset of accretion. The pic¬ 
ture emerging from the studies is that, although accretion 
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could be prompt under the right conditions, it is likely 
that the onset of accretion is delayed by an appreciable 
amount of time, perhaps of order a year. 

Our study introduces a new class of TDEs, for which 
a puffed disk or torus forms promptly after disrup¬ 
tion. Eurthermore, the BH accretes at a steady and 
highly super-Eddington rate of about lO^M 0 yr“^ for 
a few days. This h ighly super-Eddington rate resem¬ 
bles those found by Coughlin & Begelman (2014). Our 
TDEs involve ultra-close encounters between low mass 
(0.57 — 1 Mq) stars and a 10^ Mq BH. TDEs with in¬ 
termediate mass BHs, but larger separation encounters, 
have been also studied by |Ramirez-Ruiz fc Rosswo^ 
(2009). Our simulations are in a regime where account¬ 
ing tor full general relativistic effects is needed, including 
those from the spin of the BH. In these ultra-close en¬ 
counters, the star is effectively disrupted as it arrives 
at periapsis, with the tidal debris plunging into the BH 
almost instantaneously (on the timescale of one orbital 
period). Moreover, given the extreme proximity of the 
debris to the BH, general relativistic precession is very 
efficient in circularizing stellar material to form an ex¬ 
panding torus. The inner material in the torus spirals 
into the BH and is accreted at a constant rate until the 
supply of bound debris is exhausted. 


2. TIDAL DISRUPTIONS AT A GLANCE 


A star of mass M* and radius i?* approaching a BH 
of mass Mft , likely in a highly eccentric or parabolic or¬ 


bit (IRees 1988 

Magorrian & Tremaine 

1999 Hayasaki 


et aI.|l^(H2), wil 

1 be disrupted by tidal forces if it wan- 


ders within a distance to the BH 


Rt = R* 


\mJ ’ 


( 1 ) 


called the tidal radius. It is customary to characterize the 
strength of a TDE encounter by its penetration factor (3, 
which is defined as 


~ Rp Rp \M^ ) 


( 2 ) 


with Rp the periapsis distance. The fourth length scale 
in the problem is the gravitational radius Rg = G Mh/c^, 
which is equal to half the horizon radius for a non¬ 
spinning BH and the full horizon radius for a maximally 
rotating BH. 

Given i?*, Rt, Rp and Rg, it is useful to identify the 


domain of astrophysical relevance of TDEs (Luminet & 


Pichon 1989). In the (3 vs Mh plane, this domain is 


the triangle shown in Eigure[^ whose demarcations are 
obtained by interpreting (3 in Eq. as a function of Rp. 
The base of the triangle is Rp = Rf^ i.e. /3t = l3(Rt) = I- 
Below this line (Rp > Rt) there is no disruption. The 
left side of the triangle is the line obtained by setting 
Rp = R„', that is. 


I3, = (3(R,) = 



( 3 ) 


To the left of this line (Rp < i?*) lie very close encounters 
where the BH enters the star in the process of disrupting 
it. This process has many similarities to the high-speed 
collisions of stellar-mass BHs and red giants as studied by 
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Fig. 1.— Tidal disruption domain for a M* = Mq BH and 
R* = Rq star. The thick vertical line denotes ca nonical TDEs. 
The s quare point shows the ultra-close TDE from [Laguna et al.| 
||1993[l and crosses those in the present study. 


Dale et al. (2009). Finally, the right side of the triangle 


is the line when Rp = Rg; that is, 


13 a = ma) 


R* 

GM*/c2 {mJ 


( 4 ) 


To the right of this line (Rp < Rg) are events where the 
star enters the BH before it is disrupted. 

According to Figure for a solar-type star (M* = Mq 
and i?* = Rq), the maximum BH mass to disrupt the 
star is Mh = 3.2 x 10® Mq. Moreover, the maximum pen¬ 
etration factor with disruption is /3 = 78, which involves 
a Mh = 4.7 X 10® Mq BH. The edges of the triangle of 
astrophysical relevance are, of course, not sharp due to 
variations in the definition of of relevant length scales 
arising from the spin of the BH, the space-time curva¬ 
ture in the neighborhood of the BH, and the internal 
structure of the star. 

After disruption, the receding debris spreads, with 
roughly half of the material remaining bound to the BH. 
The most bound material has specific binding energy 
given by 


f-min — 


GMhAR, 


R) 


~ -G /3 C MV^Rf^Ml'^ , 


( 5 ) 


where Ai?* = ^ R* is the spread of the debris with ^ a 
deformation factor, and 

R/= y [/3n-b(I-n)] = ^/(/3,n) ( 6 ) 

is the distance to the hole when the spread in specific 
binding e nergy of the debris freeze s-in. In the original 
estimates (Evans fc Kochane^|I989 ), n = 0 (f = I), and 


thus Rf = Rr,. iVIore recent studi es (Stone et al. 2013 
Guillochon & Ramirez-Ruiz 2013) suggest that n = 


(f = (3), so Rf = Rt. In our study, we observed that 
n ~ 0.5, which for large f3 translates to / ~ 0.5. 

With emin at hand, the characteristic fallback time for 
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the most tightly bound material to return to the BH is 

GMh 

(7) 

Furthermore, from the Keplerian relation 

^ = h2TrGMhfh-^G ( 8 ) 

at 3 


and the mass per specific binding energy 
dMh M* M, 1 


de 2|e„iin| {2Tr G Mh) 


2/3 


(9) 


which is assumed to be roughly constant, the accretion 
rate is estimated to be 


Mh = 


dMh dMh de 


dt 


de dt 


^ M„ 


t 


-5/3 


^ > ( 10 ) 


with Mmax 
in Eq. (10 

= M*/(3 tmin)- The power-law decay of t 3/3 
) is considered to be a ubiquitous property 

of TDEO 

Elees||1988l lEvans & Kochanek 

1989 

Laguna 

et al.|1993| Rosswog et al.||20U9). Slight departures from 


this power-law have been found close to the peak ac- 
creti on rate as a result of the equation of state of the 
star (Lodato et al. 2009| or effects from the spin of the 


massive 


BH jHaas et al. |201 2|). 

chon & Hamirez-Kuiz ( |2ijl3[ ) has shown that the canon¬ 
ical t~^di characterizes only full tidal disruptions. 

In the case of a TDE with M* = Mq, R.^, = Rq and 
Mh = 106 Mq, 

tmin ^ 0.11 / 3"6 /3^-3/2 rl^'^ yr , ( 11 ) 


Recent work by Guillo- 


and 

- 0.3 /33 ml Mg-'/" Mq yr"', (12) 


where Me = M/i/lO^M©, r* = R^,/Rq and m* = 
M^/Mq. For comparison, the Eddington accretion rate 
in this situation is MEdd = 0-02 Mq MQyr~^ (assuming 
10 % efficiency). 


3. NEW TDE REGIME 

We are interested in TDEs with /3 = 10 and 15, involv¬ 
ing a BH with mass Mh = 10® Mq . They are denoted 
by crosses in Figure and are closer to the BH-enters- 
star boundary than the “canonical” scenarios of /3 = few 
and Mh = 10 ® Mq, which are denoted by a thick vertical 
line in Figure In the same fi gure, the square poin t 
shows the ultra-close TDE from Laguna et al. (1993). 
Here, we consider non-spinning BHs and BHs with spin 
a/Mh = ±0.65. The sign denotes whether the spin of 
the BH is aligned (plus) or anti-aligned (minus) with 
the orbital angular momentum of the star. The disrup¬ 
tions involve main-sequence stars injected in parabolic 
orbits with masses M* = 1 Mq and 0.57 Mq, modeled as 
polytropes with equation of state p oc Pg and F = 4/3. 
During the evolution, we use a gamma-law equation of 
state p = poe(r — 1). Table [^provides the parameters of 
the simulations: the penetration factor (3, the BH spin 
parameter a, and the mass of the star M*. 


The simulations were carried out with our Maya code 
and fully account for general relativistic effects. This 
code was als o used in our pre vious general relativistic 
TDE study (Haas et al. 2012). The Maya code is 4th 
order accurate in solving the Einstein equations and 2nd 
order in the hydrodynamics equations. The general con¬ 
vergence of the code is slightly above 2nd order. We em¬ 
ploy nine levels of mesh refinements. All but the coarsest 
mesh have 70^ grid-points, with the coarsest having 1203. 
The resolution on the finest mesh is Rg/50. 

In general terms, the TDEs studied here involve stars 
comparable in size to the BH, i?* ~ 4.7Mg"' r* Rg, peri- 
apsis distances of about 


Rp ~ 4.64/3fo' Mg'/^ mj'/^ R,, (13) 

and tidal radius 

~ 218Mg-"/^m/'/^r*i?g, (14) 


where Mg = Mh/10^ Mq and /3io = /3/10. Although 
the value of Rp suggests that the star could potentially 
swing by the BH without the BH entering the star, the 
combination of large (3 and general relativistic effects pro¬ 
duce an outcome dramatically different from the situa¬ 
tions with ,5 ~ 1 and 10® Mq mass BHs (next section). 

With /? > 10 and intermediate mass BHs, the star will 
be effectively disrupted and stretched to a few times its 
original size by the time it reaches periapsis passage. A 
deformation factor ^ ~ 4 was found to be c omm on in our 
simulations. Therefore, from Eqs. ([^ and (10) 




-3 j3t-3/2 3/2 -1 , ,1/2 

c ' r» m* Mg' s, 


(15) 


and 

M„,ax 9.6 X 103 Pig /-3 2 ^-3/2 ^2 ^^- 1/2 ^ 

(16) 

where ^4 = ^/4. For reference, the Eddington accretion 
rate in this case is MEdd = 0.002 Mg Mq yr”^. Our sim¬ 
ulations show that Mraax ^ 10"' Mq yr“3 and t„iin 25 s. 
This time-scale is comparable to the circular orbital pe¬ 
riod of the most bound material: 


Rc.irc. — 2 TT 


l{Rp-AR,)3 


V GMh 

~ 316 (1 - AR^/RpfG ^-3/2 ^3/2 ^- 1/2 g 

~16.5^-''/"r*/"TO)r'/"s, (17) 

where ~ 0.86 ^4/ 3 io m);/^ Mg/^. 


4. ANATOMY OF A DISRUPTION 

We focus this discussion chiefly on the accretion rates 
onto the BH. Figure!^ shows the accretion rates of tidal 
debris through the BH horizon. The time axis is such 
that t = 0 s denotes periapsis passage. Three distinct 
accretion epochs or stages are identifiable in most of the 
cases. 

The first phase is a narrow spike or flare in the accre¬ 
tion rate. The spike is due to the portion of the stellar 
debris that immediately plunges into the BH. An esti- 
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TABLE 1 

Simulation Parameters and Accretion Rates. 


Run 

P 

a/Mh 

M^/Mq 

Mmax (M0yr“l) 

Mate(A7©yr i) 

PioSqMi 

10 

0 

1 

3.6 X 10^ 

1.0 X 10® 

PioSo.esMi 

10 

0.65 

1 

1.0 X 10-* 

7.5 X 102 

/3ioS-o.65Afi 

10 

-0.65 

1 

1.2 X 10^5 

2.0 X 10^ 

PioSoMq,^'^ 

10 

0 

0.57 

3.5 X 10^ 

1.0 X 10^ 

P10S0.65M0.57 

10 

0.65 

0.57 

7.6 X 10^ 

4.0 X 10^ 

P10S-0.65M0.57 

10 

-0.65 

0.57 

7.9 X 10-* 

3.0 X 10^ 

P15S0M0.57 

15 

0 

0.57 

2.0 X 10-* 

4.0 X 10^ 

P15S0.65M0.57 

15 

0.65 

0.57 

8.8 X 10-* 

3.5 X 10^ 

P15S-0.65M0.57 

15 

-0.65 

0.57 

3.8 X 10® 

2.0 X 10^ 


mate of this accretion rate is given by 

^ -^h Poo ^oo 


' 10^ /3 


1/2 . ,7/3 
10 ^5 


7/6 -7/2 , , 

r, Mgyr 


(18) 


where we have used Ah AttRI Voo - (GMh/Rpy/^ 
and Poo ~ M*/(47ri?*/3). This estimate is consistent 
with the values for Mmax reported in Table [d 

After the accretion flare there is a decay ;^ase, which 
for a non-spinning BH (see top panel in Figure^ loosely 
resembles a power-law decay. The duration of this decay 
seems to depend on the spin of the BH and the penetra¬ 
tion factor /3. 

Finally, as new material returns to the BH, it circu¬ 
larizes and forms an accretion torus on a times cale of 
lOtmin ~ 1,400s (Evans & Kochanek 1989), with 
the accretion eventually reaching steady state at about 
Afiate 10^ Mq yr“^, as noted in the last column of Ta¬ 

ble The steady state accretion will cease when the 
mass supply from the bound debris is depleted. 

5. SIMULATION RESULTS 

To illustrate the effect of the mass of the star, panels A, 
B and C in Figure show the BH accretion rate for /3 = 
10, grouping the cases with the same BH spin parameter. 
Since the accretion rate seems to be insensitive to the 
mass of the star, we will focus on the M* = 0.57 Mq star 
given that we have a wider variety of simulations for this 
star. 

In Figure]^ we also organize the runs according to /3, 
wit /3 = 10 in panel D, and /3 = 15 in panel E. The 
peak accretion rate of the flare is higher if the BH is 
spinning. Interestingly, the case with a counter-rotating 
orbit (i.e. BH spin anti-aligned with the orbital angu¬ 
lar momentum) yields the largest peak accretion rate. 
For the /3 = 10 simulations, the post-flare accretion rate 
depends on the BH spin magnitude but not its orien¬ 
tation, which is consistent with similar findings fo r the 
steady-state, su bsonic accretion onto a moving BH (Pet- 
rich et al. (1988). The late-time accretion rate seerns to 
increase with the spin of the BH. Another difference be¬ 
tween the a/Mh = 0 and the a/Mh = ±0.65 cases is that 
for the latter, the time scale for flare decay is shorter; 
that is, the late constant accretion phase, which signals 
the formation of the torus, is reached after a couple of 
hundred seconds. 

The corresponding /? = 15 cases show that the spin 
of the BH does not play a role in determining the late¬ 
time, constant accretion rate. Furthermore, as with the 
/3 = 10 cases with spinning BHs, the constant accretion 



10'^ 

10^ 

10^ 

10^ 

IQi 


.1 


/3ioSoMo.57 

/3lO'S'-0.65Afo.57 

^lO'So.esA/o.s? 




/3l5SoAfo.57 —■ 
^155'-0.65Afo.57 ■■■• 
t*15'So.65A7o.57 -. 



1500 


Fig. 2.— BH accretion rate as a function of time. Panels A- 
C display cases with penetration factor 0 = 10 grouped by the 
spin parameter of the BH, with a/Mh = 0, —0.65 and 0.65 corre¬ 
sponding to panels A, B and C, respectively. The remaining panels 
display accretion rates for M* = 0.57 Mq with penetration factor 
/3 = 10 (panel D) and 0 = 15 (panel E). 
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Fig. 3 . — Snapshots of the density (top row), temperature (middle row) and specific entropy (bottom row) of the stellar debris for the 
/3 ioSoMo. 57 case. Columns from left to right depict times t = —28 s, 450 s and 1,250 s, respectively. 


phase is reached within a few hundred seconds. For both 
the /3 = 10 and /3 = 15 cases, the late-time accretion rate 
is ^ 10^ Mq yr“^. 

In summary, we found that the accretion rate is not 
very sensitive to M*, and that the rate reached during 
the flare depends on both a/Mh and /3. Furthermore, 
the late accretion rate for /3 = 15 is not sensitive to 
the spin of the BH, and for both penetration factors is 
approximately M ^ 10^ Mq yr“^. With the exception of 
the cases with /3 = 10 and a/Mh = 0, 0.65 (panels A and 
C in Figure]^, there are no hints of a power-law decay 
rate. 

Figurel^shows snapshots of the density (top row), tem¬ 
perature (middle row) and specific entropy (bottom row) 
of the stellar debris for the ,Sio«S'oAfo .57 case. Columns 
from left to right depict times t = —28 s, 450 s and 
1,250 s, respectively. Self-intersection of debris and sub¬ 
sequent formation of the accretion torus are evident. The 
increase in temperature and entropy as a result of shocks 
due to self-intersection of material is also evident. 


Finally, Figure shows histograms of the mass per 
unit binding energy as a function of binding energy for 
the same case and times as in Figure It is evident that 
dMh/de is not flat, as is required to achieve the 
p ower-law decay. 


Carter & Luminet (19821 were first to suggest a tidal 
compression perpendicular to the orbital plane oc 
that would yield an increase in the central density and 
temperature of the star oc and /3^, respectively. We 
do not observe such extreme compression. As mentioned 
before, the star arrives at periapsis effectively disrupted, 
with its central density decaying and undergoing only a 
mild compression, adequately handled with the resolu¬ 
tion used in our simulations. 


6. CONCLUSIONS 

We presented a new class of TDEs showing prompt 
formation of an accretion torus and hyperaccretion. 
These TDEs involve ultra-close encounters with a Mh = 
10^ Mq BH. The accretion rates are highly super- 
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Fig. 4.— Histograms of the mass per unit binding energy as a 
function of binding energy for the case in Figure]^ 


Eddington. Additionally, there is no evidence of a 
decay in the accretion rate. This is likely due the strong 
influence of general relativistic effects in these ultra-close 
encounters. The late-time accretion rate, once the torus 
has formed, reaches an approximate steady-state and re¬ 
mains highly super-Eddington at Miate ~ lO^M 0 yr“^. 
With this rate, the BH should be able to accrete the ma¬ 
jority of the bound tidal debris in just a few days. A 
word of caution is needed regarding the accretion rate 
quotes in this study. Our simulations did not include 
effects from radiation, which could potentially decrease 
the rates. However, it is not likely that radiation ef¬ 
fects would diminish the rates enough to make them sub- 
Eddington, since the Eddington rate in these situations 
is MEdd = 0.002 Ms Mq yr“^ (assuming 10 % efficiency). 

In a subsequent study, we will expand the parame¬ 
ter space of our simulations and investigate the emission 
properties of the tidal debris. It will be important to 
investigate further the role of shocks in heating and cir¬ 
cularizing the debris as well as the overall role of cooling 
effects. Moreover, we will consider inclined orbits relative 
to the spin axis of the BH in order to com pare our re¬ 
sults with those of Shen & Matzner (2014). We are also 
interested in the possibility of amplification of magnetic 
fields and whether the class of TDEs in the present study 
could provide a n explanation for jetted events, such a s 
Swift J1644-I-57 ( Burrows et al.|20il Levan et al.|2011 ). 

We thank Tarnara Bogdanovic. Work supported by 
NSF 1205864, 1212433, 1333360. Computations at 
XSEDE TG-PHY120016. 
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